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SUMMARY 

As a resxilt of a fundamental investigation of the meteorological 
conditions conducive to the formation of ice on aircraft and a study 
of the process of airfoil thermal ice prevention, previously derived 
equations for calculating the rate of heat transfer from airfoils in 
icing conditions vere verified. Knowledge of the manner in which 
water is deposited on and evaporated from the surface of a heated 
airfoil was expanded sufficiently to allow reasonably accurate calcu- 
lations of airfoil heat requirements. The reseeirch consisted of 
fll^t tests in natural— icing conditions with two 8-foot— chord heated 
airfoils of different sections. Measurements of the meteorological 
variables conducive to ice formation were made simultaneoiisly with the 
procurement of airfoil thermal data. 

It wa,8 concliided that the extent of knowledge on the neteorology 
of icing, the Impingement of water drops on airfoil siirfaces, and the 
processes of heat transfer and evaporation from a wetted airfoil 
surface has been Increased to a point where the design of heated wings 
on a fundamental, wet— air basis now can be undertaken with reasonable 
certainty. 


INTRODUCTION 

For a period of several years, the National' Advisory Committee 
for Aeronautics has conducted research on the prevention of ice 
formation on aircraft through the use of heat. DiJring this time, 
research of a fundamental nature on the problem of thermal ice 
prevention was reteirded by the more urgent need for development of 
ice— prevention systems for specific airplanes in military service. 
Satisfactory wing— and tall— evirface thermal ice- prevention systems 
for a Lockheed 12-A, Consolidated B-24, Boeing B-I7, and 
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Ciirtiss-Wrl^t C— 46 airplanes (references 1, 2, 3> and 4, respectively) 
were designed, fabricated, and tested in natural-icing conditions. 
Windshield thermal ice-prevention systems which proved adeqmte in 
the icing conditions encountered were provided for the 12-A, B— 24, 
and C— 46 airplanes. Each wing- and tall— surface design was based on 
establishing, for fll^t in clear^lr conditions, a surface- 
teD5)erature rise above free— stream ten^ieratvire which experience in 
simulated— and natiiral— icing conditions had shown to be adequate for 
ice prevention. This empirical method, while proving satisfactory 
for the airplanes tested, was limited, since it was not established 
on a fvindamental basis, and a more basic procedure foimded on 
designing for the conditions existing in icing clouds was needed. 

The MCA at present is engaged in an investigation to provide a 
fundamental understanding of the pt^ocess of thermal ice prevention 
in order (l) to establish a basis for the extrapolation of present 
limited data on heat requirements to meteorological anti fli^t condi- 
tions for which test data are not available, (2) to provide data for 
improving the efficiency of thermal ice- prevention equipment, ani^ 

(3) to provide a wetr-alr, or meteorological, basis for the preparation 
of design specifications for thermal ice— prevention equipment. The 
resesoTch consists of an investigation of the meteorological factors 
conducive to icing, and a study of the heat— transfer processes which 
govern the operation of thermal ice— prevention equipment for airfoils 
and for windshield configurations. 

The airfoil heat— transfer phase of this investigation consisted 
of the measurement of the factors affecting the transfer of heat 
from airfoil surfaces during flight in natiural— icing conditions. 

These data are correlated with the slmulteineoue measxiremente of the 
mateorologicaJ. parameters which influence the heat— transfer process, 
and are analyzed for the purpose of establishing a we1>-air ice- 
prevention design basis for airfoils. 

Qlie first approach to the icing heat— transfer problem on a 
fu n da m entetl basis was made in Englemd by Hardy and Mann prior to 
1942, In this study, a method for the calculation of heat tremsfer 
from a heated surface subjected to icing conditions was presented 
and substantiated by measurements in an icing tunnel. Later work by 
Hardy in which these heat— transfer equations were modified for 
general application is presented in references 5 and 6. Reference 5 
contains infoimation on the protection of eQ.1 aircraft conqponents 
against ice accretion. Reference 6, prepared dtaring a period of 
active participation by Mr. Hardy in the MCA icing research program, 
presents an analysis of the dissipation of heat in conditions of 
icing from a section of the heated wing of the C— 46 airplane 
(reference 4). 
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Other research in the present NACA investigation has been 
reported in references 7> 8, and 9» Eeference 7 gives the first 
measurements in this program of the liquid-^ter concentration in 
clouds# Eeferences 8 and 9 deal with the meteorologiceil aspects of 
icing conditions in stratus clouds and in precipitation areas of the 
warm-front type. 

Ee search on the problem of heat transfer from airfoils in condi- 
tions of icing has also been conducted by other laboratories. In 
reference 10, the transfer of heat from surfaces subjected to icing 
conditions on Mount Washington has been studied. The General Electric 
Eesearch Laboratory has conducted a number of investigations on this 
phase of icing. A summary of this work and a list of reports Is 
presented in reference 11. A camprehenslve report by the Army Air 
Forces on the development and application of heated wings is contained 
in reference 12. 

In continiiatlon of the present icing program, the C-h6 airplane 
was equipped with special meteorological and electrically heated test 
apparatus, and flown in natural— icing conditions diuring the winters 
of 19^5—46 and 19^^-^7» Fll^t tests were conducted mainly along 
airline routes over most of the United States. The meteorologiceLi 
data recorded during the icing conditions encountered in the two 
seasons are presented and discussed in references 13 and lU. 

This report presents an analysis of the data obtained during the 
19k5-46 and 19^6-47 winter seasons with two electrically heated air- 
foil sections. The data were analyzed tislng the heat— transfer 
equations developed by Heirdy. (See references 5 a n d 6.) A considera- 
tion of the area and rate of water impingement on one of the airfoil 
sections based on an analytical study of water-drop trajectories 
(reference 15) is also presented. An attempt is made to further the 
knowledge of the process .of airfoil thermal ice— prevention. 

The appreciation of the NACA is extended to Uilted Air Lines, 

Inc., the United States Weather Btireau, and to the Air Materiel 
Command of the Army Air Forces for aid and cooperation in the research. 
In particular, the services of Major James L. Murray of the Air 
Materiel Command, Army Air Forces, and Captain Carl M. Christenson 
and First Officer Lyle W. Eeynolds of United Air Lines, who served as 
pilots of the research airplane, were a veiluable aid to the conduct of 
the investigation. 
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SYMBOLS 

TTift following nomenclature is used throughout this report: 

radius of water drop, feet 

airfoil chord length, feet 

specific heat of air at constant pressure, Btu per pound, 
degree Fahrenheit 

specific heat of water at constant pressure, equal to 1 Btu 
per poijnd, degree FeOu^enheit 

concentration factor, defined in equation (6), dimensionless 

saturation vapor pressvire with respect to water, millimeters 
of mercitry 

water-drop collection efficiency, defined in equation (lO) 

acceleration of gravity, equal to 32.2 feet per second, second 

convective surface heat— transfer coefficient, Btu per hour, 
square foot, degree Fahrenheit 

toteuL surface heat— transfer coefficient, Btu per hoiir, square 
foot, degree Fahreniieit 

mechanical equivalent of heat, equal to 778 foo-^-pounds per 
Btu 

thermal conductivity, Btu per second, square foot, degree 
Fahrenheit per foot 

dimensionless drop-inertia quantity, defined in equation ( 5 ) 

latent heat of vaporization of water at surface teii5)erature, 
Btu per pound 

liquid— water concentration of icing cloud, pounds of water 
per cubic foot of air 

weight rate of water-drop impingement per unit of surface 
area, pounds per hour, square foot 

weight rate of water flow aft of area of watei^-drop impinge- 
ment per foot of span for one side of airfoil, pounds per 
hour, foot 
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Ru 
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Wa 

Ws 
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X 
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weight rate of watei>-drop impingement per foot of span for one 
side of airfoil, poionds per hour, foot 

concentration of liquid water contained in drops of each size 
in a drop-size distribution, pounds of water per cubic foot 
of air 

barometric pressure, millimeters of mercury 

Prandtl number (cpu/k), dimensionless 

unit rate of heat flow, Btu per hour, square foot 

Reynolds munber for airfoil (Vcy/l-i), dimensionless 

free— stream Reynolds number of water drop relative to speed of 
airfoil (2Va7/i-i), dimensionless 

distance measured chordwiso along airfoil surface from stagna- 
tion point, feet 

tenqierature , degrees Fahrenheit 

local velocity Just outside bo-undary layer, feet per second 
free— stream velocity, feet per second 

weight rate of evaporation of water per unit of surface area, 
pounds per hour, square foot 

weight rate of evaporation of water per foot of span for one 
side of airfoil, pounds per hotu’, foot 

distance measured chordwise along airfoil chord line from 
zero— percent chord point, feet 

evaporation factor, defined in equation (P2), dimensionless 

airfoil ordinate, feet 

starting distance of water drop above projected chord line of 
airfoil, feet 

pressure altitude, feet 

ratio of saturated to dry adiabatic lapse rates 


3 exponent of Prandtl number, 1/2 for laminar flow, I /3 for 

turbiilent flow 
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7 specific wel^t of air, pounds per cutlc foot 

7v specific wei^t of water, equal to 62 »k pounds per cubic foot 

viscosity of air, pounds per second, foot 
Subscripts 

1 refers to conditions at edge of boundary layer 

k kinetic 

o refers to free— stream conditions 

me mean effective 

s refers to conditions at airfoil surface 

SL sea level 


ANALYSIS 

During fli^t in icing conditions a heated wing is cooled by 
convective heat transfer, by evaporation of the water on lie surface, 
and, in the region of droplet interception, by the water striking 
the wing»^ The rate at which heat must be supplied in order to 
maintain the wing surface at a specified temperature is, therefore, a 
function of the rates of convection, evaporation, and water in^Jingement, 
Equations for expressing this heat requirement are presented in 
references 5 and 6, These equations, with slight modification, are 
used throughout this report. 

Expressed as an equation, the unit heat loss q from a partially 
or couqpletely wetted surface exposed to icing conditions may be stated; 

<1 = <lw + Q.C + 0.6 (l) 


where 



0.W 

heat 

loss 

Ic 

heat 

loss 


heat 

loss 


Each of these individual heat flows will be analyzed. 


^The heat loss due to radiation is small and can be neglected 
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Heat Loss Du© to Wanning tha Intercepted Water 

In the region where water droplets strike the wing, the heat 
required per tinlt area to heat the water to surface ten^jerature Is 


q.w = ^ 




( 2 ) 


The term Atj^^ Is the kinetic temperatxjre rise of the water caused 
hy stoppage of the droplets as they strike the wing. The value of 
Atk;, is given hy 


Atk^ 


2gJcpv, 


( 3 ) 


where V is the free— stream velocity in feet per second. The value 
of Atk^ is less than 2 ° Fahrenheit for airplane speeds up to 200 
miles an hour and, for the calculations presented in this report, the 
term has "been neglected. Equation ( 2 ) thus hecomes: 


Q.W = ^ (ts— to) 


(h) 


The weight rat© of water inqplngement on the wing, the area of 
impingement, and the distrlhutlon of the water over that area are 
important factors in the heat— transfer analysis. In addition to the 
effect of the amoimt of water intercepted on the value of q^ in 
equation ( 4 ), the ©valuation of provides an indication of the 
quantity of water which must he maintained in a liqviid state until 
it either evaporates or runs off the trailing edge if the formation 
of ice aft of the area of impingement, normally termed "runhack," is 
to he avoided. The area of impingement Influences the extent of 
heated region to he provided at the leading edge, while knowledge of 
the distrlhutlon of water iD5>lngement is required in the calculation 
of the heating requirement in areas where water is striking. 

Calcvilations have been mad© hy Glauert (reference I6) for the 
trajectories of water drops about cylinders and an airfoil. In this 
work the assumption was made that the drops obeyed Stokes* law of 
resistance. At the speeds of flight, however, Stokes* law no longer 
strictly holds, and Langmuir and Blodgett (referenj© I7) con5)ut©d a 
series of drop trajectories about cylinders, spheres, atxI ribbons, 
taking into consideration deviations from Stokes* law. Thee© con^jit- 
tations were undertaken on the assun^tion that the trajectories for 
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cylinders would apply to airfoils if the airfoil were replaced hy an 
"equivalent" cylinder (reference 12). 

Preliminary calculations "based on references l6 and 17 indicated 
that, for large values of drop size and airspeed, the assumption of 
the equivalent cylinder woifLd not hold for airfoils. Therefore, more 
extensive calculations were undertaken to determine the drop trajec- 
tories for one of the test airfoils of this research, an NACA 0012 
airfoil at 0° angle of attack. In these calculations, presented in 
detail in reference 15/ a Joukowski airfoil (the contoxrr of which 
closely approximates that of the NACA 0012) was used to supply the 
stream lines since the Joukowski stream lines and velocity field can 
he computed with relative ease. The basic eqixations presented in 
reference l6 were used with modifications for deviation from Stokes* 
law as given in reference 17 . The procedure followed was to start a 
given distance forward of the airfoil and calculate the paths of the 
drops using a step-hy— step integration process. Results of these 
computations are presented in figure 1. The curves shown establish 
the distance s, measured from the stagnation point, at which a given 
drop will strike the airfoil when starting a distance y© above 
the projected chord linb. Curves are presented for various values 
of K, where 


It should be noted that the curves of figure 1 apply strictly only 
for a drop Reynolds nixmber Ry of 95»^5/ that is, only for partic- 
ular combinations of drop size, airspeed, altitude, and air ten^jera- 
ture. The value of 95*65 was chosen as being the Reynolds number 
corresponding to average conditions of drop size, airspeed, altitude, 
and air temperature experienced during the tests of this investiga- 
tion. However, the curves of figure 1 can be used for a range of 
Reynolds numbers on either side of 95*65 without serious error. Due 
to practical considerations, these curves were used in the analysis 
of the data presented in this report, even thou^ the Reynolds 
number differed scmewhat for every case. 

Area of water impingement.— The end points of the curves shown 
in figure i denote the extreme location at which drops of a partic- 
ular K value will strike the airfoil. Beyond this value of s/c no 
drops of this K value will hit. Thus, the broken line in figure 1 
establishes the area of impingement for all values of K. 

Rate of water impingement .- The rate of water impingement at a 
specified point on an airfoil is a function of the area of impingement, 
the velocity of flight, the liquid— water concentration of the air 



(5) 
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stream, «rii^ the dlstrlhutlon of the intercepted water over the surface, 
This latter factor, called the concentration factor C is represented 
hy the ratio of y^ to s, or: 


For point values. 



or more exactly. 


The wel^t rate of water ln5)ingement per unit of surface area in 
pounds per hour, square foot, then, is 


It is apparent from equation (6) that C is simply the slope of 
the curves shown in figure 1. A plot of the measured slopes of these 
curves as a function of s/c is presented in figure 2, Using values 
of C obtained from figure 2, the wei^t rate of water iiq>ingement 
at any point on the surface cem he calculated from equation (7)» 

In the case of a cloud, where the water drops are not of uniform 
size, hut instead follow a pattern of size dlstrlhutlon, the late of 
loplngement can he confuted if the dlstrlhutlon is known or assumed. 
Th e rate of water impingement at any point is the' sum of all the rates 
of Impingement of the volume of water contained in each drop size, 
Eqimitlon (7) then hecomes 


where n is the concentration of liquid water cont ain ed in drops of 
a partlciQar size C Is the concentration factor for the K 

value coriTespondlng to that drop size. 

In order to estahllsh the posslhlllty of runhack forming aft of 
the heated area of a wing, it is necessary to know the total quantity 
of water Intercepted per unit of wing spcm. This rate of impingement, 
denoted as Mb in pounds per hour, foot span, is given hy 



( 6 ) 


a 3600 T m C 


(7) 



(8) 
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Mb 



ds 


( 9 ) 


A mare rapid msthod for the evaluation of Mg utilizes a curve of 
collection efficiency E as a function of K (fig. 3 ). Collection 
efficiency is defined as 


E - yoiimlt 

ymai 


( 10 ) 


where yoiimit value of jo for which drops of a particular 

K veilue Just miss the airfoil, and ymax is the max imu m ordinate of 
the airfoil. The equation for can^iuting Mg, then, is 


Mg = 3^00 EVm Virwr (ll) 

Using figure 3, the rate of water impingement can he conputed for 
each of the drop sizes in the assumed or measured drop-size dlstri— 
hutlon. The total rate of impingement is the summation of these 
individual rates. 


Heat Loss Due to Forced Convection 

The unit heat flow from the surface of a body in an air stream 
resulting from convective heat transfer ceui he expressed: 


qc = h (tg— tojc) 


( 12 ) 


where tg is the surface tenq)erature emd tot is the kinetic temper- 
ature of the free-etream air at the point for which the heat flow is 
being computed, ihe factor h is the convective heat-transfer 
coefficient and may he evaluated by measurements in clear air or by 
calculation using the methods presented in references l8 and 19 * 
Evaluation of the term toj^ will now be discussed. 

The sxirface of an unheated wing moving throu^ the air will 
assume a temperature somewhat hi^er than that of the free air stream 
because of stoppage of the air particles in the boimdary layer next 
to the surface. This temperature rise is of importance in the calcxj- 
lation of heat requirements for ice prevention in that it establishes 
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the datum point frcaa >diich the teniperature of the sva*face must he 
raised to obtain the desired temperature, ts. The value of the 
temperature rise in clear air, from equations derived in reference 
is, for laminar flow. 


( 13 ) 


and for turbulent flow. 


where U is the local velocity Just outside the bou n da r y layer at the 
point along the sxjrface where the value of Atj^ is being calculated* 

In clouds, the kinetic temperature rise is induced, due to 
evaporation of water fl*om the surface. Assuming the surface is 
completely wetted with water, the value of the temperature rise for 
laminar flow becomes 



Atk 


2gJcp 


[i - ^ (i - 


( 15 ) 


where 


©1 



(16) 


and eok is the vapor pressure at saturation at the wet kinetic tem- 
perature, tok* The value of tok is 


or 


tok = + Atk 
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= to + 


Jt- 

2gJ‘cp 


_ IE 

Y2 


1 \ 1 
- Pr^ J 


- 


\~^J 


( 17 ) 


Equation (I7) Is for ]ftmlna.-r flow. The equation for turtulent flow 
Is the same, hut with the exponent of Pr changed to l/3» It can he 
seen that this equation must he solved hy trieil, since the value of 
eojj. is dependent upon the temperature toj^. 


Experiments in clouds, in the process of calibrating a free-air 
thermometer installation (reference I3), showed that hy nrultiplying 
the clear-air kinetic— temperature rise hy the ratio of the saturated 
to the dry eidlahatic lapse rates, good agreement between the values 
of kinetic temperature rise calculated in this manner the measured 
values was obtained. Since use of the ratio of the adiabatic lapse 
rates was substantiated e35>erimentally, and since eqviation (17) must 
be solved by trial, a somewhat laborious procedure, the following 
equations were used in this report to calculate values of toj^: 

For laminar flow. 






(18) 


and for turbiilent flow. 


tojj. 


to 


+ JC2- 

2gJcp 





(19) 


Values of ae/ctd> the ratio of the wet— to the dry— adiabatic lapse 
rates, are obtained from figure 4. The use of the lapse— rate ratio 
in equations (I8) and (I9) is seml-empirlcal. The limitations of 
this slnq>lification in the calculation of kinetic— temperature rise of 
airfoil surfaces in clovids are not known. Below speeds of 200 miles 
per hour, however, these equations can be used with small error, since 
the kinetic— temperature rise is low. 


Heat Loss Due to Evaporation of the Water on the Surface 

The amount of heat removed from a wetted surface as a result of 
the evaporation of water on that surface can be expressed: 

9e “ ^s ^a 


( 20 ) 
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From reference 6 the relation "between and the conTectlve 
heat— transfer coefficient h can "be expressed for a can 5 )letely 
wetted surface as: 


where 


Qe = h (X - 1) (tg-toj^) 


X 


1 + 




Cp 


\te— ti 




Ofe 


0.622 

Pi 


( 21 ) 


( 22 ) 


By substituting average valties for I»s cp, equation (22) can 

"be rewritten 


^ ^ 3.75 (23) 

The values chosen for Lg and cp are 1100 Btu per pound and 0«24 
Btu per pound, ^Fahrenheit, respectively# The factor Psl/Pi Is "the 
ratio of the standard sear-level pressure to the local static pressure. 

It should "be noted that the evaporation factor X applies only 
when the surface is con 5 >letely wetted. If only partial wetness 
prevails, the value of X must "be modified according to the degree 
of wetness. 


Total Heat Loss from a Wetted Surface 

Summarizing the heat losses due to water in 5 )lngement, convection, 
and evaporation, equation (1) can "be written: 


q ~ ^ ("tg— to) + h (tg— to]^) + h (x — 1) (tg— tojf) 


which reduces to 


q = 1% (tg— to) + h X (tg— tojf) 


( 24 ) 
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Aft of the region of water lB5)lngeB»nt, * 0 and eqxiatlon (24) 
hecames 


q_ = h X (ts— to^) (25) 


DESCRIPTION OF EQUIPMENT 

All tests reported herein were made in the C-46 airplane shown 
in figure 5« ^Eh© airplane had "been modified to provide thermal ice- 
prevention equipment for wings, en5>ennage, windshield, and propellers. 
A description of the thermal system for the wings enq>ennage is 
given in reference 20, The windshield system was altered for the 
flights as descrlhed in reference 21, Protection for the propellers 
was provided hy electrically heated "blade shoes. 

Ihe meteorological equipment tised dxiring the tests to measure 
the free— air temperature, llq\ild-water concentration, drop size 
drop— size distribution is descrl"bed in references I3 and 14, 

Two electrically heated test airfoils were used to obtain 
fimdamental data on the process of wing thermal ice prevention. Each 
airfoil was mounted vertically on top of the fuselage of the C-46 
airplane, as shown in figure 5* The test airfoil Installed during 
the winter of 1945—46 had an NACA 0012 section. For the tests in the 
winter of 1946-47 the airfoil had an NACA 65,2— 016 section in order 
to provide test data for low-drag sections, as well as conventional 
sections. Both sections are symmetrical, cuad the models were 
Insteilled with the chord line in the plane of symmetry of the air- 
plane; that is, at zero single of attack for \uiyawed flight. Ordinates 
for an NACA 0012 airfoil are given in reference 22, and for an 
NACA 65,2— 016 airfoil, in reference 23, Figure 6(a) shows the 
NACA 0012 airfoil mounted on the fuselage, the NACA 65,2—016 
airfoil was mounted as shown in figure 6(b). A clear plastic 
blister, shown in figures 6(a) and 6(b), allowed the airfoils to 
be viewed and photographed in flight. 

Both airfoils had an 8— foot chord and a 4,7— foot span, with a 
faired square tip, A heated test section of 1— foot span was located 
2 feet above the top of the fuselage. It had been determined 
previously, by means of a pressure survey, that the test— section 
location was well above the edge of the fuselage boundary layer, 
ElectrlcaiJLy heated guard sections were built auround the leading- 
edge region on both sides of the test section for the purpose of 
preventing any disturbaince of the air flowing over the test section 
which might have been caused by ice accretions in the region of the 
guard sections. 


MCA TN No. 1472 


15 


NACA 65,2-016 Airfoil Model 

Construction details of the 65,2— OI6 airfoil model sire shown in 
figure 7. The metal portion of the structirre consisted of aluminum 
rihs and skin supported from the fuselage hy two spars. The test 
section was made up of a 3/8-inchr-thick plastic "base and a sheet of 
plastic— in5>regnated fabric, l/6It— inch-thick, on top of which 1/2— inch 
wide, 0,002— inch-thick, electrical resisteince heating strips were 
cemented in a spanwise direction spaced I/32 inch apart, A covering 
of the l/6I»— inch-thick plastlc-ii!5)regnated sheet was laid over the 
resistance strips, and on top of this was cemented a skin of 
0,006-lnchr-thick aluminum. Each l/ 2 — inch— wide heating strip was 
connected to individual lugs located along the edges of the test 
section. This provided means for chordwise adjustment of the power 
distribution by l/ 2 — inch increments. The heated area of the test 
section extended back to 77 percent chord on the left side and to 
17 percent chord on the right side. 

The guard sections were constructed in the earns manner as the 
test section, with the exception that the aJ-uminum skin was 0.011 
inch thick. The heated area of the guard sections extended to I7 
percent chord on both sides of the model, 

Msasurements of the temperatxn^ of the alianlnum surface of the 
test section were obtained by means of thermocouples. Fine iron- 
constantan thermocouple wire was rolled flat to produce a strip 
approximately 0,002 inch thick and I/16 inch wide. These strip 
thermocouples were laid in spanwise grooves about 3 inches long cut in 
the aluminum skin. The thermocouple Junctions were located in the 
middle of the grooves, and the leads passed through holes at the ends 
of the grooves into the interior of the model. Aluminum was sprayed 
into the groove over the strip thermocouple for a distance of about 
3/16 inch on either side of the Junction, Thus, the thermocouple 
Junction was bonded to the alxmiinum skin, allowing accurate surfaoe— 
temperature measuremsnts to be made. The remainder of the groove on 
either side of the aluminum spray was filled with a nonelectrically 
conducting material. Thermocouples were located at the center of 
the test section at 1— inch chordwise intervale in the leading-edge 
and ceilculated transition regions, and at 1—1/2— inch chordwise 
intervals in other regions. Surface temperatures were recorded by 
means of self- 4 )alancing automatic— recording potentiometers. 

The flow of heat through the outer surface was calculated from 
measiirements of the power dissipated in the electrical heating 
strips. This power was determined by measuring the resistance of 
the strips the cxurrent flowing through them. Thermocouples 
placed on both surfaces of the plastic base at a number of chordwise 
stations gave an indication of the heat flow into the ^del interior. 
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These thermocouples were connected to the same recording potenti- 
ometers used to record siarface ten^jeratxires. 

Pressure taps were installed flush in the test-section surface 
about 3 inches down from the top edge of the test section at varioxis 
chordwise points for the purpose of neasuring sxcrface pressure dis- 
tribution. A standard NACA 60 -cell pressure recorder was xised to 
record the pressures. 

A source of 400 -cycle, single— phase, alternating current was 
supplied to the test and guard sections for heating these sxnrfaces. 
The heating strips for the test section were grouped into 3O chord- 
wise blocks. Control of the current flowing throu^ each block was 
provided, so that a large variation in the chordwise distribution of 
heat flow was possible diuring flight. Before the icing operations 
steirted the heating strips for the guard sections were connected to 
maintain a constant sxirface— ten5>erature rise during flight in clear 
air. Controls were provided so that the total heat input to the guard 
sections, but not the chordwise distribution, could be varied during 
flight. 


NACA 0012 Airfoil Model 

With a few exceptions, the construction of the 0012 airfoil 
model was substantially the same as that of the 65,2— 016 model. These 
exceptions will be noted. 

The top layer of plastic-in^jregnated fabric covering the electri- 
cal resistance strips constituted the outer skin of the test and guard 
sections. This was painted and sanded, after the test-section thermo- 
couples had been installed. The heated area of the test section 
extended back to 58 percent chord on the left side and to 11 percent 
chord on the right side. The heated area of the guard sections 
extended to 11 percent chord on both sides of the airfoil. 

Strip thermocouples of the same type as installed in the 65,2-016 
model were used to measure sxirface ten5)eratures of the test section. 
Spanwise grooves were cut in the plastic-in^jregnated fabric sheet at 
various intervale along the chord. The strip theimocouples were laid 
in the grooves and cemented in place. The surface was then painted 
and sanded so that only a thin layer of paint covered the thermo- 
couple Junctions. Thermocouples were located at the center of the test 
section at 1-inch chordwise intervals in the leading-edge and calcular- 
ted transition regions, and at 2— to 2— 1/2— inch chordwise intervals in 
other regions. Surface ten5)erature8 were recorded by means of an 
automatic— recording potentiometer. 

Die flow of heat throu^^ the surface of the test section was ceJ.- 
culated from measurements of the power dissipated in the electrical 
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haatlng strips, in a manner similar to that descrlhed for the 65,2—016 
airfoil. 

Installation of surface pressure taps for the measurement of 
p3ressure distribution was the same as for the 65,2—016 model. The 
pressvcres were recorded by photographing a multiple— tube manomater 
to which the pressure taps were connected. 

A source of 400-cycle, single— phase, alternating current was 
supplied to the test and guard sections for heating these surfaces. 

In the test section, provisions were made for obtaining a limited 
mamber of chordwlse heating distributions as well as for control of 
the total heat input with each distribution. A small degree of variai- 
tion of each heat distribution was also provided. During fll^t it 
was possible to control only the total heat input, and to vary, to a 
small extent, each distribution. As with the 65,2— 0l6 airfoil, the 
heating strips for the guard sections were connected to give an 
approximately consteint surface-temperature rise in clear air. No 
control of the heat distribution or the toteJ. power input to the guard 
sections was provided. 


OEST PROCEDURE 

The test airplane was flown into natinral- icing conditions over 
most of the northwestern area of the IMited States during the winter 
of 1945—46. Ihrrlng the winter of 1946-47 the area of operations was 
extended to Include a few flights in the central euad eastern part of 
the Iftilted States. The usual test procedure, during flight in icing 
conditions, was to record> airfoil data simultaneously with the 
measurement of the meteorological conditions. The rotating cylinders, 
described ip reference I3, which constituted the means of measuring 
llquld-^water concentration, drop size, and drop-size distribution, 
were extended as often as was conveniently possible. Records of 
free-air tengjerature, airspeed, and altitude werb taken several times 
a minute. The recording potentiometer vised to obtain airfoil temperei— 
tvu:es was operated continuously. During this time, the values of 
current flow through the electriccQ. heating strips of the airfoil 
were recorded. Photograph? of the test-section surface and records 
of pressure distribution were taken at frequent intervals. 


RESULTS 

A tabulation of the fll^t and meteorological conditions for 
which simultaneous airfoil data were obtained is presented in tables I 
and II. Table I contains the fll^t and icing conditions for which 
corresponding heat— transfer measurements were made with the NACA 0012 
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airfoil, Tatle II gives similar Information for the NACA 65,2— 0l6 
airfoil. All measurements were made dttring flight In natureil— Icing 
conditions. During most of the flights, large variations in liquid- 
water concentration, and occasionally drop size, were experienced. 

The rotating cylinders, used to msasure liquid— water content and drop 
size, were extended for about 1 to 2 minutes, thus giving average 
values for 1 — to 2 -minute intervals. A complete cycle of the 
recording potentiometers used to record airfoil tenqjeratures required 
h to 6 minutes, dxirlng which time the meteorological conditions may 
have changed considerably. For these reasons, an effort was made to 
select, for analysis and discussion herein, only the airfoil data 
recorded during flight in relatively uniform clouds euid/or where 
close correlation existed between the cylinder meas\ireiBenta and the 
airfoil— temperature records. Of these data, only a part, chosen as 
being typical, are presented in this report. These are the thermal 
data for which the flight and icing conditions are given in tables I 
and II. 


NACA 0012 Airfoil Data 

Figures 8(a) to 8(g), inclusive, present the measurements of 
surface temperature, surface heat flow, and resulting heat— tiansfer 
coefficients obtained with the 0012 airfoil model during flight in 
the conditions presented in table I. The heat— flow distribution 
Illustrated in these figures had been found by experiment to give 
an approximately imiform temperature rise over the test-section 
surface during fli^t in clear air. Variations in the intensity of 
the distribution for the different conditions of table I occurred 
as a result of the heat supply procedure followed during the tests. 
In general, dxrring an icing test the total heat input was reduced 
until the siirface temperature was observed to fall close to freezing 
temperature at some point on the test section, lypiceil values of 
surface temperatm'e, heat flow, and convective heat— transfer coef- 
ficient obtained during flight in clear air are shown in figure 8(h). 

The data presented in all figures except figure 8(g) were taken 
with the entire test section heated. Figure 8(g) presents data 
secured with only the leading-edge region heated, from 11 percent 
chord on the rlg^t side to 8 percent chord on the left side. At the 
time of this test, insufficient heat was supplied to the leading- 
edge area to evaporate all the water striking the surface, and 
streamers of ice formed aft of the heated region, similar to those 
shown in figure 9» 

The heat— flow values given in figures 8(a) to 8(h), inclusive, 
were calciilated from measurements of the total power dissipated in 
the electrical resistance strips and the Internal heat loss. The 
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measurements of surface ten5)erature, which were obtained with the 
thermocouple Install I at Ion previously described, were corrected for 
errors incxared by the presence of a layer of paint covering the 
function. Magnitude of jbhese errors was determined from knowledge 
of the thickness eind thermal conductivity of the paint and the amoxant 
of heat flowing through it. 

The kinetic tenq^rature of the free-stream air toj^ used in 
computing the heat— transfer coefficients, was ced.culated for the 0012 
airfoil from eq^uatlons (l8) and (19) > using experimentally determined 
values of the expression 


1 - ^ (1-PrP) 

where 3 is l/2 for laminar flow and I/3 for turbulent flow. Values 
of this expression for various points along the airfoil surface, were 
obtained from figure 10, which presents data obtained during flight 
in cleeur air. 

A typical record of pressure distribution over the 0012 airfoil 
‘model test— section svirface is shown in figure 11, 

KACA 65,2-016 Airfoil Data 

Figures 12(a) to 12 (j), inclusive, present the measurements of 
surface tenqperature, sxn^ace heat flow, and resulting coefficients 
of heat transfer obtained with the 65,2-016 airfoil model during 
flight in the conditions presented in table II. The distribution 
of heat flow shown in these figures had been experimentally estab- 
lished to provide an approximately uniform tenqperature rise above 
free-air teinperature over the test— section surface during flight in 
clear air at an altitude of 11,000 feet and a true airspeed of I75 
miles per hour. Ihls heat— flow profile was used throu^out all the 
fll^t tests. Slight variations in heating Intensity are due to 
variations in internal heat flow and chardwlse heat condxictlon in 
the thin aluminum skin. Typical values of surface temperatoore, 
heat flow, and convective heat— transfer coefficient for fllj^it in 
clear air are given in figure 12 (k). 

The results shown in all figures (except figs. 12(h) to 12(j), 
inclusive) were obtained with the test section heated to approximately 
55 percent chord. Heating aft of this point was precluded by a 
malfunctioning of the heating equipment in this area. Figures 12(h) 
to 12(j), inclusive, present data secured with only the leading-edge 
region heated, from I7 percent chord on the right side to I7 percent 
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chord on the left side. During the time vhen the measurements of 
figures 12(i) and 12(j) were taken, an insufficient quantity of heat 
was "being supplied to the leading-edge region to evaporate all the 
water striking the surface, and ice accumulated aft of the heated 
area. This is shown in figure I3. 

The heat— flow values given in figures 12(a) to 12(k), inclusive, 
were calculated from measurements of the total power dissipated in 
the electrical resistance strips, in a similar manner to that used 
for esta"blishing the heat flow for the NACA 0012 airfoil. In 
addition to the determination of the internal heat loss in computing 
the surface heat flow for the KACA 65,2— 016 airfoil, the flow of 
heat chordwlse in the thin alumlnxan surface was considered. The 
chordwise heat conduction is a fxmction of the chordwlse variation 
in surface temperature. It was assianed that the qxiantity of heat 
flowing from point to point along the surface, as indicated "by the 
difference in surface temperatiure between the two points, originated 
from the heating strip under the hl^er temperature and flowed away 
from the surface into the air stream in the area of strface over the 
heating strip at the lower temperature. This method, although 
Inexact, offered a rapid means of estimating the effect of chordwlse 
conduction. A more exact determination of this effect can be , 

obtained using the "relaxation” method of reference 10. No correc- 
tions were applied to the surface— temperature measurements, since 
it was assumed that the surface thermocouple Junctions were at 
sxirface tenqieratxire . 

The kinetic tenq>erature of the free— stream air toj^ was 
calculated vising equations (I8) and (19)» "Values of the expression 
1 — (1— were calculated and are plotted in figvure l4. 

A picture of the conditions of wetness >diich existed on the 
airfoil during flight in clouds cem be seen in figure 15. This 
figvire shows some typical records obtained with strips of blueprint 
paper which had been fastened to a device that could "be extended 
into the air stream up the dealing edge of the airfoil model to a 
point Just below the test section. Since, in effect, these were 
wrapped arovind the leading edge of the model, they illustrate the 
pattern that the water assumes in striking the airfoil and flowing 
aft. The records were obtained during icing conditions 11, I3, and 
14, table II. 

A typical measurenent of pressure distribution over the 65,2—016 
airfoil model test— section surface is shown in figure I6. 
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DISCUSSION 

The ultimata In performance of a wing thennal ice— prevention 
system is one idilch will prevent the accretion of ice on any portion 
of the wing. This ideal operation requires that any water on the 
wing sxn^ace must he maintained in a liquid state xmtll it evaporatee, 
or blows off the wing at the trailing edge. In many wing designs^ 
heating of the entire surface is not practicable because of such 
features as Integral fuel tanks^ arxl in these Instances any water 
flowing aft of the heated region is apt to freeze and form the type 
of Ice accretlcai normally termed runback. In the following dis- 
cussion the ice— prevention action of a heated wing will bo examined 
in detail, and the reliability of the equations arwi assimtptlons 
presented in the analysis section for the prediction of surface 
tenqeratxire and rate of water evaporation from the surface will be 
established. If these equations and analytical methods can be shown 
to define correctly the process of thermal ice prevention, the funda- 
mental design procedure for a heated wing initially conceived in 
reference 6 will be more firmly established. The enq)lrlcal design 
method of providing a specified temperature rise in clear air can 
then be replaced by the more fundamental smd flexible concept of 
supplying sufficient heat to maintain the surface temperature above 
freezing \mtll the water is either evaporated or carried away. 

An analysis of the action of a heated wing requires the consid- 
eration of three factors: namely (l) the meteorological and fll^t 
conditions for which the wing must provide protection; (2) the area 
of water l2i5)lngement, and the rate and distribution of impingement 
over that area; and (3) the rate at which the water is evaporated 
from the wing surface. 


Meteorological and Flight Conditions 

The specification of a meteorological condition for the design 
of thermal ice-prevention equipment depends upon the geographical 
areas over ^ich the airplane will fly, the seasons of opeivitlon, and 
other factors dictated by the Intended service of the aii^raft. 
Obviously, the establishment of design conditions for a specific area 
requires a knowledge of the conditions prevailing over the area. If, 
on the other hand, the ice— prevention system is to provide protection 
for all-weather operation, general specifications of a meteorological 
condition must be established which will encoiopass all conditions 
likely to be encountered. 

The most recent and extensive information in regard to the 
severity of icing conditions likely to be e:^rienced in all-weather 
operation in the IMited States is contained in references I3 and l 4 . 


22 


NACA TN No. Ihl2 


In reference 13, estimates of the maximum continuous icing condi- 
tions as well as the maximum probable icing conditions apt to be 
encountered are presented. Since the duration of the maximum 
probable icing condition is quite short (1 to 2 minutes), and icing 
of this severity is entirely associated with cumulus clouds which 
should be avoided in all operations, the maximum continuous icing 
condition is believed to be of greater interest for design purposes. 
Two conditions of maximum continuous icing are presented based on a 
relationship of drop size and liquid-water content. These conditions 
are given in the following table: 


Liquid-water 

concentration 

(gm/m®) 

^Msaik-effective 
drop diameter 
(microns) 

Free-air 

tenq)orature 

{°F) 

0.8 

15 

20 

0.5 

25 

20 


It is believed that the conditions in the above table form a good 
basis for the design of thermal ice— prevention equipment for eill— 
weather operation. In addition to these values, however, the proposed 
wing thermal system shoxild be analyzed for possible undesirable oper- 
ation in other icing conditions. For exaiiq)le, reference I 3 points out 
that drops of 35 to 50 microns diameter should not be regarded as 
exceptional. Although the amount of liquid water associated with such 
large drops is usu a ll y low (about 0.1 gram per cubic meter) the fact 
remains that tlie area of water 'impingement would be very large 
would probably exceed the limits of the heated region if this region 
had been based only on a consideration of the data in the maTiTnnm 
continuous table. Finally, the possibility of encountering icing con- 
ditions at low tenq>eraturee may be a critical condition for heated 
wings on some ali^lanes. For insteuice, the estimated conditions 
of maxlmaam continuous icing presented in reference I 3 given 
in the preceding table were extended in reference Ih to air tenqera- 
tures as low as —20® F for the case of 15-mlcron drops. The con- 
ditions of maximum continuous icing suggested are 0.5 gram per cubic 
meter at 0° F and 0.25 gram per cubic meter at -20® F, both with 
2 — — 

The msan-effective diameter as defined in reference I 3 is the size of 

drop in a cloud sanq>le for which the amount of liquid water existing 
in water drops larger than that drop is equal to the amoaant of liquid 
water existing in drops smaller them the drop. 
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a maan-effectlve drop diameter of 15 microns. Either of these cond.i— 
tions may he more d.eleterious to the functioning of the ving thermal 
system than those in the table. 

From the foregoing discussion, it is evident that the analysis 
of a heated wing should give consideration to several possibly criti- 
cal icing conditions in the same manner that several flight condi- 
tions are assumed for the ving structural analysis. The data of 
references 13 and l4, although somewhat limited in scope, are consldr- 
ered to be ST;ifficlently indicative of icing conditions in the Nbited 
States to form a meteorological basis for heated-wing design. 

The problem of selecting a fli^t condition for the design of 
ice— prevention equipment is concerned with the airspeed and eLLtltude 
at \dilch the airplane will fly. The airspeed will depend igpon the 
specific airplane, and, in general, a cruise condition should be 
selected. Choosing an altitude for design is dependent upon several 
factors, \diich will be discussed later. 


Area, Bate, and Distribution of Water Iiiqplngement 

Having defined the icing conditions for idilch the heated ving is 
to be designed, the next step is to determine the region of the 
leading edge in which the water drops will strike the wing, the rate 
of water impact at any specified point in that region, and the total 
rate of iiiq>lngem9nt per foot of wing span. This subject was dis- 
cussed at some length for the genered case in the analysis section. 

In that discussion, it was shown that the method of reference 15 could 
be used to prepare (for any wing section for which the stream lines 
were known or could be determined) cturves similar to those presented 
in figures 1, 2, and 3. Ohe broken line of figure 1 gives an Indicar- 
tlon of the area of water in^lngement, while the rate of tnq>lngement 
at a specified point can be obtained from figure 2 and eqmtloin (7)» 

Two methods are available for the determination of the total 
rate of impingement per foot of wing span. The calculation of this 
q^uantlty is of primary lDq>ortance, as It detezmlnes the amount of 
heat required to disperse the water by evaporation. The first method 
utilizes the concept of collection efficiency E as mentioned In 
the analysis section. Oils method Is preferable idien only the value 
of the total rate of impingement Is desired, since preparation of 
the curves of figure 2 Is not required. For a thorou^ analysis of 
the heat transfer from the surface, however, knowledge of the rate 
of lBq>lngement at a point Ma Is required. By enqtloylng eqtiation (7) 
and figure 2 a curve of the distribution of water In^lngement 
(Ma against s/c) ccui be plotted. Figure I 7 shows such a curve for 
the NACA 0012 airfoil, using equation (8) and an **E" type drop-size 
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distribution. (See reference I7.) A cvarre of this type presents an 
Interesting plctxire of the distribution of water iaqplngement and, in 
addition, the area under the curve denotes the total rate of water 
interception. 

Althou^ the method of reference 15 is considered to provide a 
con^jlete and q.ulte accurate prediction of the distribution of water 
la^dngement on the leading edge of an airfoil, it does have the dis- 
advantages of requiring (1) a knowledge of the velocity cofl 5 )onents 
along a number of the airfoil stream lines, and (2) considerable 
conqjutatlon. The difficulties associated with the con^iutatlon of 
the water-drop trajectories for airfoils have encouraged the substi- 
tution of a cylinder with reidlus equal to the airfoil leading-edge 
radlTis in the determination of water ioq>lngement. (See references 5 
and 6.) The curves of reference I 7 , which have been calculated for 
a large range of drop sizes, airspeeds, altitudes, and cylinder 
diameters, are then \ised directly to evaluate the anticipated water 
lB 5 )ingement on the airfoil. This substitution proced\n*e is a useful 
device but should be employed with a full knowledge of its limita- 
tions. One of these limitations is the fact that the curves of 
reference I 7 provide the area and total rate of water impingement, 
but give no direct indication of the distribution of tn 5 )ingement, 

A second restriction of the cylinder-substitution method is 
concerned with the contour and size of the forward portion of the 
airfoil. To obtain an indication of this effect, the rate and area of 
water ln5)i]3gement on the 0012 airfoil, at 0° angle of attack, and on 
the leading-edge cylinder of that airfoil are conq)ared for the same 
flight conditions and various drop sizes in figure I 8 . The values for 
the 0012 section were obtained frcan figures 1 and 2, and those for the 
cylinder from reference 17. At drop diameters up to about 25 microns the 
rates of impingement on the airfoil and on the leading-edge cylinder 
are approximately the same, eilthou^^ above 25 microns as the drop size 
increases, the rate of irqjlngement becomee considerably greater on 
the airfoil than on the cylinder. At drop diameters up to about I 8 
microns the area of Impingement on the cylinder is roughly equal to 
that on the airfoil. However, at a drop diameter of 25 microns, -(diich 
is not unusual (reference I 3 ), and was presented previously in this 
report as a possible wwtI rmim continuous condition, the area of inplnge— 
lasnt on the airfoil is nearly 50 percent greater than on the cylinder. 

It should be noted that the value of 25 microns for the maximum contin- 
uous condition is the mean-effective diameter, and that drops of a 
larger size probably will be present dtie to the existence of a distri- 
bution of sizes. Althou^ these values provide an indication of the 
scale limitation of the cylinder— substitution method, the fact shoxild 
be noted that figure I 8 applies to only one airfoil section, with 
an 8-foot chord, and at one flight condition. The leading-edge radius 
of the NACA 0012 section for an 8-foot chord is small (1.5 in.) 
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aod the leadlng-ed^ cyllader does not seitch the section contour for 
any great extent shore the chord line. This Is shown graphically In 
flgrnre 19 which presents a con 5 )arlson of the forwaard portions of 
three airfoil sections and the leading-edge cylinder of the 0012 
section. In the case of airfoil sections with the leeidlng-edge radliis 
a greater percent of the chord than the 0012 section, anil also for 
airfoils of 0012 section, or similar, with chords greater than 8 feet, 
the cyllnder-suhstltutlon method will present a better approximation 
than that Indicated hy figure l8 for the same speed range. 

for airfoil sections with a leading edge radius idilch represents 
a Bma l l percentage of the chord, the substitution of an ”equlTalont" 
cylinder (ref essence 12) with a reidlus larger than the leading-edge 
radltis would probably provide a better Indication of the rate 
area of water In^lngement on the airfoil than would be obtained for 
the leading-edge cylinder. At the present time there Is not suffic- 
ient Information on watei'-drOp trajectories about airfoils to provide 
a basis for selecting the proper cylinder In each Instance; therefore, 
the designer must utilize the more coiQ>llcated, but moz*e acctirate, 
method of reference 1^ or assume some cylinder diameter based on his 
ejqperlence. The possibility that the rate ar><^ area of Infringement on 
an ellipse would more closely approximate the rate and area of 
Impingement on a series of similar airfoils has been siiggested and 
Is worthy of future cooslderatlon. 

The ability to select a proper drop size for the design of wing 
Ice-^reventlon equipment Is a fcu:tor of considerable Importance to 
the designer, as can be Illustrated by figure 20. In this figure the 
rate and area of impingement are presented for the 0012 airfoil as 
a function of drop size. The rate of Infplngement for each drop size 
was calculated for a liquid-water content of 1.0 gram per cubic meter. 
Consider, then, a change In design drop diameter from 10 microns to 
20 microns. Ttie resultant Increase In rate of water laglngement Is 
1.75 pounds per hour per foot of span or an IncroMe of I 75 percent, 
although the actual amount of water present per \mlt yoltaae of cloud 
has not been changed at all. The same Increase In drc^ size will 
cause an Increase in area of lo^lngement from 1.^ to 4 percent s/c. 

In contrast, consider the effect on the rate of water ladings— 
sent produced by an Increase In the quantity of liquid water present, 
assuBiing the drop size to remain constant. The area of lsq)lngement 
will remain imchanged, idille the rate of lsq)lngement will Increase 
only In direct proportion to the increase In water concentration. 

This exasqple clearly illustrates the fact that the asiotmt of free 
water p^sent In an Icing cloud Is only one factor influencing the 
quantity of water which will actually strike the wing In a specified 
tlsm Interval, that the size of the cloud drops is a factor of 
at least eqiial Isqportance. 
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The prohlem of dlstrihution of the sizes of drops In an icing 
cloiid also hears careful consideration. For example, if an "E” type 
drop-size distribution exists with a nean-effectJ-^ drop diameter of 
25 microns, the largest drops will he 68 microns In diameter. Hence, 
the area of impingement will he conslderahly greater than if a 
uniform drop size of 25 microns prevailed. The data of reference l4 
indicate that, in general, the distrihutions of drop size in icing 
clouds are fairly narrow, ai3d do not usvially follow the hroad dlstri— 
hutlons, such as type "E.” Nevertheless, the distrlhution must he 
considered, since the largest drops in the cloud determine the area 
of lD 5 )lngement and the minimum extent of heated area required for 
ice prevention. 


Rate of Evaporation of Water 

Having discussed the problem of area and rate of water intercep- 
tion, the next step is to establish the rate at which the Intercepted 
water is evaporated from the wing surface emd the validity of the 
equations presented previously for determining the rate of heat 
dissipation during the process of evaporation. The problem of rate 
of evaporation is particularly lm 5 )ortant because all of the water 
intercepted by a wing heated only in the region of the leading edge 
must be dispersed by evaporation if the formation of runback is to 
be avoided. 

From a superficial study of the mechanism by which water is 
deposited on the surface of a wing, it woxild be expected that in the 
area of water lnq)ingement the surface is completely wetted, and that 
equation (2k) for calculating the heat loss from a heated wing is 
valid. of the area of impingement, it would be anticipated 

that the surface may not be fully wetted, since water does not reach 
this region directly, but instead must flow back from the area of 
impingement. If the surface aft of the area of ii^ingement is only 
partieilly wetted, the expression for X (equation (23)) must be 
modified for tise in equation (25) to calcxilate heat requirements. 

Observations, made d\u:lng the airfoil tests, of the water 
pattern on the leading edge of the airfoils revealed that the above 
suppositions are correct. At a very short distance back of the 
region of ln^iingement, the film of water was observed to reach a 
state of instability and break into small rivulets. A picture of 
the conditions of wetness which actually exist on a wing dtirlng 
flight in clouds ceui be seen in figure I 5 . This figure, which shows 
typiceO. records obtained with the strips of blueprint paper placed 
the leeuilng edge of the 65,2— OI6 airfoil model during flight 
In icing conditions. Illustrates the pattern formed by the water in 
strikii^g the airfoil leading edge and flowing aft. It is evident 
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from the patterns that the area of ln^ilngement, vhich is clearly 
defined. Is con5)letely wetted, while hack of this area the water 
collects and forms rivulets, creating a partially wetted surface. A 
study of the patterns indicated a variation in the fraction of surface 
area wetted (aft of the impingement area) with rate of impingement of 
water. Accordingly, the rates of water la5>lngenient Mg were calcTi— 
lated using eq^uation (11) for the conditions existing at the time 
that rivulet patterns were obtained (icing conditions 11 through 15, 
table II). The ciirves and values presented previously for the 0012 
airfoil were used in the calculations of Mg for the 65,2— 0l6 airfoil. 
Substitution of the calculations for the 0012 section in computing 
values for the 65,2—016 section appears to be a good approximation, 
since the contour of the 65,2—016 section in the leading-edge region 
is very nearly the same as that of the Joukowskl airfoil used in the 
0012 trajectory calculations. Figure I9 coiapeLres the contours of 
the three sections. 

The values of Mg were plotted against the measured areas of 
surface wetted, obtained from the strips of blueprint paper. 

Figvire 21 shows the relationship, thus obtained, between the rate of 
flow of water from the impingement region and the fraction of stirface 
area wetted. For the data shown in figure 21, values of the rate of 
water flow over the surface aft of Implngemsnt Mb were assumed 
to be equal to the rate of water impingement Mg. Ihe scatter of data 
points in the figure is believed to be caxxsed by errors in measure- 
ment of the liquid— water concentration occitrring at the time the 
rivulet patterns were obtained. Table n shows that the free-air 
temperature was high during icing conditions 11 throu^ I5, when the 
blueprint records were taken. The kinetic tenperatvire was close to 
freezing, gmd it was observed that the water striking the rotating 
cylinders, used in the measiirement of water concentration, was 
running back, eind possibly off, the cylinders. Thus, the liquid- 
water concentrations measured may have been lower th^ the actual 
concentrations present. The two data points corresponding to a 
wel^t rate of water flow of 0.57 pounds per hour per foot of span 
(fig. 21) represent the rivulet data procured at the lowest free- 
air tenperatvire of these tests (icing condition lU, table II). These 
are probably the most reliable data, since the rotating cylinders 
were subject to smaller losses of water. Hierefore, the curve shown 
in figure 21 was weighted toward these points. The ultimate extent 
of this curve in the direction of percent of svirface wetted is not 
definitely known. Ihere is evidence, however, indicating that the 
degree of svirface wetness eift of the area of impingement reaches a 
maximum which is not exceeded, regardless of the rate at which water 
is intercepted. It was found that a relationship exists between the 
rate of flow of water in the region eift of impingement and the 
svirface— temperature rise above free— air temperature, and that the 
temperature rise decreases to a limit as the rate of flow of water 
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Increases. Figure 22 shows the relationship "between the rate at 
which water flows "back over the heated surface of the 65,2— OI6 aii>- 
foil test section and the average Increase in ten5>erattare a"bove 
free-air temperature of the surface from 10 to 25 percent chord. The 
values of the rate of water flow were o"btained "by subtracting from 
the calcxilated rates of water liiQ>ingement for one side of the airfoil 
the com^ruted rates of evaporation from the region of in5)ingement. 

Figure 22 iUxxstrates that as the rate of water flowing over the 
surface increases, the ten^eratvore of the surface decreases, hut that 
a limit to the decrease in teuqperature apparently is reached. This 
indicates that the rate of e-vaporation reaches a maximum as the 
limiting surface ten^jeratTure is approached, since evaporation is 
the only variable in the heat— transfer process in the area of surface 
under consideration. Therefore, if the rate of evaporation attains 
a maximum, the degree of surface wetness must also approach a limit. 

It can be demonstrated, using equation (25) and the values presented 
in figures 21 and 22, that the maximum fraction of surface area 
wetted is about 50 percent. 

Althoii^ the data from which the curve of figure 21 was computed 
were obtained with blueprint paper strips wrapped around the leading 
edge of the 65,2— 016 airfoil, the values given in this figure are 
believed to be sufficiently Indicative of the conditions of wettability 
existing on ell clean wing surfaces not specially treated to be appli- 
cable for general airfoil thermal design. For purposes of design, 
it is suggested that the limit of surface wetness for surfaces not 
specleilly treated be taken as 40 percent. It is of in5>ortance to 
note that in tusing the curve of surface wetness shown In figure 21, 
for a heated wing, the total rate of evaporation of water Wg in 
the region of water-drop impingement must be subtracted from the 
total rate of water inqplngemsnt Mg in order to obtain the rate of 
flow of water rearward from the area of la^jingement. The values 
given in figure 21 for degree of surface wetness are believed to be 
accurate only to the nearest 10 percent. 

With the information gained so far, it should be possible to 
analyze the data obtained with the two electrically heated airfoil 
models and establish the validity of equations (24) and (25) for 
calculating heat flow. The ciurves of measured heat— flow distribution 
shown in figures 8(a) to 8(g) and figures 12(a) to 12(j) were faired 
to produce a form more suitable for comparison with heat— flow curves 
calcvilated using equations (24) and (25). Conq)arlsons of the 
neasured heat flow and the heat flow calculated to produce the 
measured surface temperatures, assuming the entire surface to be 
completely wetted, for the two airfoil sections for typical cases are 
shown in figures 23 and 24. These curves are also compared to the 
calculated heat loss due to convection only, that is, assuming the 
surface to be completely dry (equation (12)). In the previously 
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mentioned calctilatlons^ measured ralues of convectlTo heat— transfer 
coefficient, obtained during flight In clear air, were used. In 
order to calculate the amount of heat dissipated In warming the 
Impinging water (equation (4)) for the calculation of heat flow for 
the con^iletely wetted case, walues of were computed from equa- 

tion (8) using the measured values of liquid— water concentration, 
drop size, and drop-size distribution (tables I and IX). As was 
done previously in the analysis of the blueprint— paper rivulet 
patterns, values of were confuted for both airfoil sections 
using the cxirves presented for the 0012 airfoil. 

A study of the measured and calculated heat— flow ctnrves In 
figures 23 and 24 shows that In the area of wate]>-drop lBq>lngement 
good agreement Is obtained between measured values and the values 
calculated for a completely wetted surface. Indicating that In the 
region ^ere It Is reasonable to assinne a fully wetted surface the 
equations for calculating heat flow are valid. Aft of the area of 
Impingement, In the region of low heat flow, iriiere it has been shown 
that the surface is only partially wetted, the values calculated for 
a completely wetted siirface are lower than the measured valties. Since 
the sin*face is only partially wetted, it would be expected that the 
calculated curve, idilch represents the values of heat flow required 
to produce the measured surface temperatures If the surface were 
completely wetted, would be considerably higgler than the mieasured 
curve. There appear to be only two possible explanations for this 
discrepancy: (1) equation (2^) gives erroneous values and cannot be 
relied upon for calculation, and (2) the valties of convective heat- 
transfer coefficient used In equation (25) for calculating the heat- 
flow values are In error. The first ejq)liaiatlQn does not appear to 
be lUcely In view of the fact that the eqimitlon was derived on a 
sound basis. Also, there Is no obvious reason why the equation 
should hold In the leading -edge region and fall to hold In the area 
aft of the leading edge. Ihe second explanation, that erroneous 
values of convective heat-transfer coefficient were used In the cal- 
culations, seemm entirely possible. Since the values of convective 
heat— transfer coefficient used In equation ( 25 ) are those measured 
during flight In clear air. It seems reasonable to assume that tran- 
sition from laminar to turbulent flow mioved forward during fll^t in 
Icing conditions from the position maintained In clear air. Movement 
of transition to a point near the leadli^g edge would cause the convec- 
tive heat— transfer coefficients In the region under consideration to 
be Increased several tlmies above the values existing In clear air, 
since the convective coefficients in tiarbulent flow generally are 
considerably greater than those In laminar flow. Svich an Increase 
In the convective coefficients would raise the curves calcxilated for 
a coopletely wotted surface (flg^. 23 and 24) to a position above the 
measured, curves. / 
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It Is very ^likely that disturljance of the hovmdary layer, caxised 
hy water drops striking the airfoil surface and roxighenlng the surface 
as they ooeOasce and flow aft, would effect a forward moveiosnt of 
transition* There is further evidence to support the assuniptlon that 
the wateiwou^ened sxrrface caused movement of transition forward* 
Observations of the 65,2-016 airfoil during a flight in clear air 
with the test section heated indicated that transition had shifted 
forward by a considerable amount* This was noted by a lowering of 
the surface ton5>eratures in the region aft of the leading edge* The 
heat distribution had been set previously to produce a constant 
surface ten5)erature in clear air, and only a change in the boimdary— 
layer characteristics could cause the evident change in heat-transfer 
coefficient* After the flight, a close examination of the leading- 
edge region of the airfoil revealed small insects stuck to the 
surface where they had hit during the flight* The surface was wiped 
clean and during a subsequent fll^t in clear air it was noted that 
transition had moved back again, an evidenced by the restoration of 
the surface temperatures to normal* Thus, it appears that very aniai i 
irregularities in the s\irface, such as are present on the sirrface of 
an airfoil in icing conditions, are sufficient to cavise transition to 
occur prematurely* Tests in wind tunnels also have shown that 
protuberances in the leading-edge region of an airfoil will cause 
the movement of tremsltion forward* (See reference 2k,) 

Most of the curves of heat— transfer coefficient measured during 
fli^t in icing conditions, shown in figures 8(a) to 8(f) and 12(a) 
to 12(g), display a definite increase in the aft region of hl^ heat 
Intensity, suggesting that transition is located at this point* It 
should be noted that the increase in heat— transfer coefficient 
indicated by these c\nrves is believed to be only an apparent Increase, 
caused by the rapid change in heating intensity in this region* If 
the coefficient is relatively constant throughout this area, as it is 
believed to be, a sudden increase in heating intensity will not be 
accoD5)anied by an equally rapid change in the thermal boundary layer, 
and for a short distance aft the indicated values of heat-transfer 
coefficient will be erroneously high* 

The exact values of the convective heat— transfer coefficient in 
the region aft of the area of inq)ingement in icing conditions are 
unknown, but it is believed the values fluctuate due to changes in 
the location of transition during fli^t* Very probably, the 
disturbance to the boundary layer caused by water on the airfoil 
surface is of such a character as to create instability in the 
boundary layer, and cause the location of transition to fluctuate* 

In the aft region of high heat flow (figs* 23 and 2k), the 
values of convective heat— transfer coefficient are known, since 
turbiilent flow existed in this region in clear air, when the values 
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were measured, as well as in icing conditions. In this region, the 
measured heat-flow curve aind the caJ-culated curve of convective heat 
traoisfer come together, indicating that at the point where the cvirves 
coincide all the water on the surface has been evaporated. 

Since the eqiaations for calculating heat flow Imve been shown to 
be valid in the area where the s\jrface is completely wetted, it is 
reasonable to assume that the equations hold in regions where the 
surface is only partly wetted, provided the correct modifications are 
made to the evaporative factor X, Fairly accurate modifications to 
the factor X are believed to be possible by using the citrve of 
starface wetness shown in figure 21. By the use of this curve it 
should be possible to cailculate the rate of evaporation of water from 
the surface aft of the region of water impingement. In the region of 
in^ingement the calculation of rate of evaporation is straiglatforward, 
since full evaporation occurs. If the rates of evaporation from the 
two test airfoils can be demonstrated to be equal to the rates of 
water inq)ingement for the test conditions, the msthod for calcula- 
ting rate of evaporation will be substantiated. 

Accordingly, calculations of the rates of evaporation from the 
surfaces of the two test airfoils were made for all the conditions 
of tables I and II for which thermal data were obtained. The rates 
of evaporation were determined graphically, using the curves of 
measured heat flow and calculated convective heat loss similar to 
those shown in figures 23 and 2k, Aft of the area of Jjiroingement, 
the position of the convective curve was established by dividing the 
measvired values of heat flow by the modified values of X (equa- 
tion (25)). Values of the degree of surface wetness used in modify- 
ing X were determined from figure 21, using con^juted rates of 
water lnroingement and evaporation from the area of impingement. The 
position of the re— calculated ciarves of convective heat loss are 
shown typically in figures 23 and 2h, The total rate of evaporation, 
then, was determined by measuring the area between the mea,sured and 
re— calculated convective heat— flow curves. This gave the total 
amount of heat dissipated by evaporation of the water in Btu per 
hour per foot of span. Dividing this value by Lg the latent heat 
of vaporization, the total rate of evaporation Wg in pounds per 
hour per foot of span was obtained. 

The rates of evaporation, obtained in the previously mentioned 
manner, are compared with the rates of water impingement, calcxilated 
by the method previously presented, for the 0012 and 65,2— 016 air- 
foil models, for the left side only, in tables III and IV. An 
average agreement of I 3 percent for all the conditions analyzed 
where no rvinbaclc formed was obtained, indicating the degree of 
reliability of the method for calculating the rate of evaporation 
of water from a heated wing. 
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In order to demonstrate further the dependahillty of the method 
for calculating rate of evaporation from a heated wing, the photograph 
of runhack on the 65,2— 0l6 airfoil (fig. 13 ) was analyzed. If it 
can he shown that the act\ial rate of formation of runhack compares 
closely to the rate at which runhack is calculated to form xmder the 
particular icing conditions, the method for calculating rate of evap- 
oration will he further suhstantiated. 

The nuiback shown in figure I3 had stai*ted forming 10 minutes 
earlier. At the time of the photograph the fonnatlon was estimated, 
hy observation during flight, to he approximately 3/1^ inch thick. 

The eirea of the formation extended about 2—1/2 inches chordwlse ari<^ 

12 inches spanwise, making a weight of ice of 0,2 pound. This 
constitutes an actueO. rate of formation of runhack of 1,2 poxmds i>er 
hour per foot span. During this 10-^nute period, two sets of 
rotating-cylinder and airfoil heat— transfer data were taken. These 
correspond to icing conditions 9 and 10, table II. Results of calcu- 
lations of the rates of impingement and evaporation based on these 
data are given in table 17, For icing condition 9 the rate of water 
in5)ingement was 1,60 pounds per hour per foot span. The rate of 
evaporation from the heated area was 0,44 poimd per hour per foot 
span, leaving a calculated rate of formation of runhack of l,l6 povuids 
per ho\ir per foot of span. During icing condition 10 the calculated 
rate of tn5)lngement was 1,79 pounds per hour per foot and the rate of 
evaporation was 0,51 pound per hour i>er foot, resiilting in a rate of 
formation of runhack of 1,28 pounds per hour per foot of span. The 
calculated rates of formation of runhack (I.16 and 1,28 Ih per hr, ft) 
agree remarkably well with the actual rate of formation (1.2 lb per 
hr, ft), illvistrating the reliability of the procedure for calcii- 
latlng rate of evaporation, 

A short time prior to this test, the airfoil was subjected to 
a much less severe icing condition (condition 8, table II), during 
>dilch all of the water intercepted was calculated to have been 
evaporated (icing condition 8, table IV). Photographs of the test 
section verified the fact that no runhack had formed. 

The foregoing analyses were based on the assimiption that removal 
of the water striking the airfoil surface is effected by evaporation 
only, and that none of the water is dispersed by mechaniceil means. 

This is consistent with the resvilts reported in reference 10 and 
25, It is believed that "blow^ff" of water, as suggested in 
references 5 and- 6, does not occur. Also, it is believed there was 
no "bounce— off" of the water drops striking the airfoil surfaces, as 
proposed in reference 11. At speeds hi^er than those encompassed 
by the scope of this investigation, it is conceivable that mechanical 
removal of the water by bounce— off co\ild occur. However, in view of 
the lack of information on this phenonenon, and since neglecting the 
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possl'bllity that water may he removed hy mechanical means tends to he 
more conservative in the thermal design, it is stiggested that bounce— 
off he neglected in the design of wiiag thermal ice-prevention equip- 
ment. 


Calculation of Heat Requirements for 
an NACA 0012 Airfoil 

Since it has heen demonstrated that the rate of evaporation of 
water from a heated wing can he calculated with reasonable certainty, 
the rate of heat flow req^llred to produce a particuleo:’ rate of evapora- 
tion can he determined with equal dependability, provided the coeffi- 
cients of convective heat transfer are known. Using the equations 
and method presented for calculating the rate of evaporation of water 
from a heated airfoil surface, a ceilculatlon was made to establish 
the extent of heated area required for ice prevention in specified 
conditions of icing for the HACA 0012 airfoil, assuming a particular 
heat— flow distribution. The conditions of calcxilation are as follows; 


Chord length ...... ..... 8 ft 

Pressvjre altitude 12,000 ft 

True airspeed 170 U 5 )h 

Free-air ten^ierature 20° F 

/ v3 

Mean— effective drop diameter ............. 25 microns 

Drop-size distribution ..E 


The procedure enq)loyed was to assume a reasonable Intensity eind dls— 
•tributlon of total heat flow and then calculate the extent of heated 
area required to evaporate all of the intercepted water. The method 
of solution will be outlined briefly in the following paragraph. A 
detailed step-by-step consideration of the problem showing all compn- 
tations is given in the appendix. 

First, the area, rate, and distribution of water impingement on 
the airfoil were calculated for the assumed conditions. Using the 
assumed distribution of total heat flow, the heat loss due to convec- 
tion for the particular conditions was then ceJ-culated. Since in an 
icing cloud the presence of water on an airfoil surface causes pre- 
mature transition, for these calculations, trsuisltlon was assumed to 
start at 5 percent s/c, and the estimated form of turbxilent heat- 
transfer coefficient shown in figure 25 was used. Calculations were 
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made for a mmiber of chordwlse stations, and the resxilts are given 
in figure 26, which shows the assumed heat— flow distribution and the 
calculated convective heat loss for one aide of the airfoil. The 
rate of evaporation is represented by the area between the curves of 
convective heat loss and total heat flow, except in the region of 
water iB5)ingement, where the rate of evaporation is denoted by the 
area between the curve of heat loss due to warming the intercepted 
water and the curve of heat loss due to convection. These areas 
actually give the rate of heat lose due to evaporation; however, by 
dividing the area value by Lg, the latent heat of vaporization, the 
rate of evaporation is obtained. The procedure, then, was to extend 
the total £ind convective curves until the total rate of evaporation 
equaled the rate of water ia^jingement. Extension of the heated area 
to l8 percent s/c was found to be adequate to ensure evaporation of 
1 of the water intercepted. 

Several other calculations were made for the 0012 airfoil to 
determine the effects of altitude, air temperature, and location of 
transition on the requirements of heat flow and extent of heated 
area necessary to evaporate all the intercepted water. The results 
of each of these calculations were compared with the results of the 
calcxilatlons for the conditions previously specified. For each of 
the calcxilations, the same total heat^flow distribution was assumed 
eind the extent of heated area required to evaporate all the intei>- 
cepted water was calculated for each condition. 

To determine the effect of eiltltude on the heat requirement, a 
comparative calculation was made for sear-level conditions with eill 
other fll^t conditions as previously specified and with the aurea 
and rate of water lBq>ingement the same as at the 12,000-foot condi- 
tion, The results of this calculation are shown in figure 27, which 
compares the relative convective heat losses at sea level and 12,000 
feet. For the conditions at 12,000 feet extension of the heated 
area to l8 percent s/c was shown previo\isly to be adequate to 
ensure evaporation of all of the water Intercepted, At sea level it 
would be necessary to extend the heated area to 26 percent s/c for 
evaporation of all the water intercepted. The curves of figure 27 
can also be used to determine the amount of Increase necessary in the 
total heat flow if all the water is to be evaporated in an area 
forward of a specified chord point. For exan5>le, assume that the 
extent of heated region for the cturves of figure 27 is limited to 
l8 percent s/c. At 12,000 feet all of the water would be evapo- 
rated, as has been previously mentioned. At sea level, however, 
some of the water would not have been evaporated. By measurement of 
the areas of figure 27 it can be shown that the total heat flow 
required to evaporate all the water within the area frcan 0 to l8 
percent s/c at sea level is approximately 10 percent greater t ha n 
the amount required at 12,000 feet. The increase in heat requirement 
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with decrease in altitude is due to the fact that the rate of evapo- 
ration of water decreases as altitude is decreased, "because of the 
decrease in the evaporative factor X (equation (23)). Since the 
convective heat— transfer coefficient increases with decrease in alti- 
tude, due to the increase in air density, it mi^t "be expected that 
the rate of evaporation would "be increased with decrease in altitude, 
"because the rate of evaporation is directly proportional to the con- 
vective coefficient (equation (2 l))» However, the increase in the 
rate of evaporation is more them conqjensated by the increase in coi>- 
vective heat loss, and the rate of evaporation, for a fixed total 
heat flow, actually becomes less with decrease in altitude. Appar- 
ently, then, airfoil thermal ice— prevention equipment in which the 
heat flow is fixed, such as electrical, systems, shoiild be designed 
for the minimum altitude at which the airplane is expected to 
encounter icing. However, if the airplane is designed to utilize 
some form of air-4ieated system, the performance of which probably 
will decrease with increase in altitude, the maximum altitude at 
which icing is expected to be encountered should also be investi- 
gated. 

To determine the effect of air temperature on the heat require- 
ment, a calculation was made of the convective heat loss at 0 F 
free— air temperature and is conq)ared in figure 27 with the convective 
heat lose at 20° F, In the calculation with the free-air temperature 
at 0° F, it was determined that the surface temperature dropped to 
freezing at 2h percent s/c before all the water on the surface was 
evaporated. However, the total heat flow required to evaporate all 
the water within the area from 0 to I8 percent s/c with the air 
temperatiire at 0° F is approximately only 15 percent greater than the 
amount required at 20° F. Althou^ this is an appreciable increase 
in the heat requirement, it is considers-bly less than that necessary 
for a similar change in conditions for ice— prevention equipment 
designed on the basis of maintaining the surface tenq)erature just 
above freezing, such as for the case of windshields. (See refers 
ence 21.) It appears, then, that a wing thermal system which has 
been designed for a relatively high air ten5>eratvire will be capable 
of ice prevention at low air tenq)eratures in icing conditions nearly 
as severe as those upon which the design was based. Of course, the 
system is more subject to failiare through the possibility of the 
surface temperature falling below freezing in the low air— temperature 
conditions, but in general, the surface temperatures required for 
evaporation of all inqjinging water in the relatively small heated 
area of the leading edge will be sufficiently high to obviate this 
possibility. 

To establish the effect of the location of transition on the heat 
requirement, a calculation was made of the convective heat loss, 
assuming laminar flow exists throughout the heated area. For this 
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calcidatloii, the neasizred values of convective heat— transfer coeffi- 
cient shown in figure 25 were used. The convective heat loss for 
laminar flow is con5)ared in figure 27 with the convective heat loss, 
assuming transition started at 5 percent s/c. In the case of con^lete 
laminar flow, it would he necesssiry to heat only to lU percent s/c 
to obtain evaporation of all the water. The total heat flow required 
to evaporate all the water within the area from 0 to l4 percent s/c 
with transition at 5 percent s/c is approximately 10 percent greater 
than the amount of heat required if laminer flow prevails. Apparently, 
the location of transition moves forward in conditions of icing, even 
in the presence of a favorable pressure gradient, to a point >diere a 
strong favorable presstire gradient is encountered (figs, 11 and l6). 

As was stated previously, the location of transition is believed to 
fluctuate, probably over a considerable disteince. It is suggested 
that foiTiard movement of transition to a point close to the leading 
edge of the wing be assumed in the design of thermal ice— prevention 
equipment, especially in view of the fact that a greater amount of 
lieat is required for the turbulent— flow condition, 

Prcm a conprehensive study of the resiilts shown in figure 27, 
some genereuL conclusions can be reached. It is apparent that aft of 
the area of droplet impingement, the efficiency of removal of water 
by evaporation decreases rapidly. The reason for the decrease in 
efficiency is that only partial wetness prevails aft of the area of 
inplngement, while the area of inplngement is entirely wet. This 
indicates that the larger the portion of the total amovuit of water 
Intercepted that is evaporated in the area of interception, the 
greater the efficiency of the thermal system becomes. The rate of 
evaporation of water is the determining factor in the efficiency of 
a wing thermal ice— prevention system. Only the heat that is dissi- 
pated in evaporation is used to advantage. The heat lost by convec- 
tion only wanns the air, Hius, the conclusion is drawn that the 
heating should be concentrated as much as possible In the leading 
edge of a wing, in the area of drop inplngement, if an efficient 
thermal system is to be obtained. 


Calculations for the C-h6 Wing Thermal System 
in Maximum Continuous Icing Conditions 

An analysis of the C-fi-6 airplane wing thermal ice— prevention 
system for the upper surface at wing station 157 was made in an 
effort to determine whether the thermal system could cope with the 
maxi.Tmnn continuous icing conditions given previously at the beginning 
of this discussion. The assumed icing and flight conditions for the 
calculations are as follows: 
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Condition A 

Condition B 

Altitude (ft) 

6000 

6000 

True airspeed (iiq)h) 

180 

180 

Liquid water 
content (gm/m®) 

0.8 

0.5 

Mean-effective 
drop size (microns) 

15 

25 

Free— air 
teB5>erature (°F) 

20 

20 


The heat— flow distrihution at station 157 estinated, "based on 
data presented in references 4 and 26, and is shown in figxare 28. 

The rates of water ia5>ingeiaent for the two icing conditions assuming 
a ”C" type drop— size distri"bution were calculated for the leading- 
edge cylinder of the airfoil section using the data presented in 
reference I7. Cu2?ves of distrihution of water impingement, for the 
upper surface, for the two cases given in figure 29, These were 
constructed using the data frcm reference 17 for each drop size in 
the distributions. The value of Mgj obtained from equation ( 9 )> 
for Condition A is O.65 pound per hour per foot spein, while the veLLue 
of Mb for Condition B is 1.39 pounds per hour per foot span. As in 
the calculations presented previously, the rate of evaporation of 
water from the surf*ace was determined by calculating the heat loss 
due to convection for the two conditions. Yalues of convective heat- 
transfer coefficient were taken from figure 30# which shows the values 
measured during fli^t in clear air using an electriceilly heated glove 
and the estimated convective coefficients for icing conditions, when 
transition moves forward. The estimated values were used in the cal- 
culations. The conqjuted curve of convective heat lose for Condition 
A is shown in figure 28 . Results of the calculations of rate of 
evaporation for the two conditions indicated that sufficient heat 
was supplied to the upper wing s\n*face to evaporate all of the water 
intercepted. For Condition A the wing was shown to be capable of 
evaporating O.9O pound per hour per foot span, indicating that the 
liquid— water concentration could attain a value of at least 1.1 grams 
per cubic meter at a mean-effective drop size of I5 microns "before 
nmback would form. The rate of evaporation for Condition B was 
calculated to be 1.35 povinds per hour per foot span, s\iggesting that 
the liquid-water concentration of O.5 gram per cubic meter is the 
limiting condition at 25 microns. 
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A further analyels was made of the upper wing surface at station 
157 using the values of Condition A, excepting that a free-air ten^ra- 
ture of 0° F was assumed. The curve of convective heat loss for this 
condition is shown in figure 28. Ihider this condition calculations 
revealed the wing would he able to evaporate O.7O pound per hour per 
foot span. The rate of impingement, as before, was O.65 pound per 
hour per foot span. 

These calciilations substantiate the general observations of the 
successful operation of the C— 46 wing thermal ice— prevention system. 

The absence of runback on the wing upper surfaces during a great many 
of the icing flights Indicates the adequacy of the thermal design. 


Selection of Conditions for Design 

In selecting values of drop size, liquid-water concentration, 
air temperature, and altitude for the design of thermal ice— prevention 
equipBBnt, a combination of these variables normally occurring in 
nature should be chosen such as to require the hipest rate of heating. 
As stated previously, conditions of maximum continuous icing are 
believed to form a good basis for design. It is of interest to inves- 
tigate the effect of different possible combinations of the variables 
of drop size, liquid— water content, air tenq^erature, and altitude on 
the heat requiremsnt for ice prevention for the maximum continuous 
conditions given in the table in the first part of this discussion. 

The effect of an increase in the size of drops in an icing condi- 
tion is to increase the collection efficiency of the airfoil, thereby 
increasing both the rate at which water is Intercepted and the area 
of impingement. An increase in the liquid— water concentration of 
the air causes a proportional increase in the amoimt of water inter- 
cepted, for a given drop size. Since all of the water striking the 
wing must be evaporated to avoid the formation of ice, an Increase 
in the rate of water Interception will cause ein Increase in the heat 
requirenent. Fortixnately, a relation between water concentration 
wnii drop size appears to exist in icing clouds, and the existence of 
very large drops generally is accompanied by a small concentration 
of liquid water (reference I3). The selection of a combination of 
drop size end water concentration should be such as to produce the 
hipest rate of impingement. It was shown previously that an increase 
in the drop size produces a greater Increase in the rate of water 
interception than a proportional increase in the water concentration. 
For this reason, the maximum continuous icing condition of 25 microns, 
rnean-effective drop size, and 0.5 gram per cubic meter, liquid-water 
concentration, generally will resiilt in a more rapid rate of water 
in^ilngemsnt than the condition of 15 microns and 0.8 gram per cubic 
meter. 
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A decrease in free-air ten^ierature, while increasing the heat 
requirement for thermal ice prevention, is accompanied hy a decrease 
in the liquid— water concentration (reference l4), which causes a 
proportional decrease in the rate of water inq)ingement, for the same 
drop size. The sizes of drops existing at low air tenq>eratures (0° F) 
In icing conditions tend to he only slightly smaller than those at 
higher air tenq^eratures (reference 13 )j therefore, a selection of 
air temperature for design will he determined hy the comhination of 
air temperature emd water concentration (and corresponding drop size) 
to produce the highest heat requirement. It will he foxmd, generally, 
that the rate of heating required to evaporate the larger qiiantlties 
of water at the higher air temperatures is greater than the heat 
needed for ice prevention at the lower tenqperatures. However, low 
aii^temperature conditions should he investigated to ascertain that 
the temperature of the heated surface will not fall below freezing. 

There appears to he no relation between altitude and the drop 
size CT liquid-water coiwentration of icing conditions, (See refer- 
ence 13 ,) Therefore, the altitude at which the heat requirement is 
greatest should he chosen. The minimum altitude of operation was 
shown previously to produce the hig^st heat requirement for wing 
thermal ice prevention. However, as was formerly suggested, if the 
aiinplane is designed to utilize some form of air-iieated system, the 
wiftTiimitn altitude at which icing is expected to be encountered should 
also he Investigated, 


C0NCLT3SI0NS 

From the foregoing discussion, it is concluded that the extent 
of knowledge on the meteorology of icing, the impingement of water 
drops on airfoil surfaces, and the processes of heat transfer and 
evaporation from a wetted airfoil surface has been increased to a 
point where the desigpi of heated wings on a fu ndame ntal, wet— air 
basis now can he undertaken with reasonable certainty. In addition 
to this general conclusion, the following conclvislons are drawn, 
based on test data and analytical studies of the processes of heat 
transfer evaporation from a heated wing: 

1, Thft heat should he concentrated as much as possible in the 
leadingj-edge region of the wing in the area of water-drop impiiigse— 
ment, if an efficient thermal system is to he obtained, 

2, An Increase in altitude, for the same rate and area of water 
inqplngpmsnt on a wing and for the same conditions of tmie airspeed 
and free— air tenqperature, decreases the heat requirement for thennal 
ice prevention. 
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3 . A wing thermal ice— prevention system which has heen designed 
to evaporate all impinging water in the leading— edge region for a 
relatively high free— air temperature (20° F) will he capable of ice 
prevention at low air temperatures (0° F) in icing conditions nearly 
as severe as those upon which the design was based. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 


APPENDIX 

Calculation of Extent of Heated Area Required 
for NACA 0012 Airfoil 

The detailed, step— by— step calcxilations for establishing the 
extent of heated area required for ice prevention on an NACA 0012 
airfoil in specified conditions of icing are presented in this appen- 
dix. It is believed that the general procedure outlined herein will 
be applicable for the design of most wing thermal ice— prevention 
equipment. 

The calculations were made for one side only of the airfoil. 

The assumed flight and meteorological conditions used in the calcur- 


lations are as follows; 

Pressure altitude . 12,000 ft 

True airspeed ....... .. 1?0 n^h 

Free— air temperature, ......... 20° F 

Liquid— water concentration. ....... 0,5 gm/m° 

Mean-effective drop diameter, ..............25 microns 

Drop— size distribution E 


The chord length of the airfoil was taken as 8 feet. 

Step 1.— Calculate area, rate, and distribution of water inter- 
ception, The area of interception is determined by the largest 
droplets present in the cloud. For the case of an "E" type drop- 
size distribution (reference 17)> "the largest drops will be 2,71 
times larger than the meein— effective size, or 
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Maximum drop diameter = 2,71 x 25 = 68 microns 

The K— value corresponding to a drop diameter of 68 microns for the 
assumed flight conditions vas calculated using equation (5). Thus, 


K = 2 

9 




^0. 


62.1;^ 


053' 


^ 4 x^. 28x10 
8 


^ 250>^X0.05^ ^ 
3 1.16xL0~® ''' 


0.465 


From figure 3^ "the efficiency of in5)ingement E for this K— value is 
54 percent. Using equation (10), the starting ordinate of the 68- 
micron-diameter drop which Just hits the surface is 


^°limit " ® 
or 

yoilmlt _ ymax 

c c 


Since the airfoil is 12 percent thick. 


^ = 0.06 

c 


and 

= 0.54 X 0.06 = 0.032 

c 

Using the broken curve in figure 1, the area of Impingement was fo\md 
to be 10.8 percent s/c. 

The distribution of water— impingement rate over the airfoil 
surface was calculated using equation (8), The indlvldml rates of 
Impingement for each of the drop sizes in the assvoned distribution 
were calciolated for various points along the surface. Ihis was 
accomplished by computing the value of K for each of the drop sizes 
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in the distribution using equation (5), determining the values of C, 
for the corresponding K— value, at various points along the svirface, 
then evaluating the expression VnC. By adding these individual 
irapingement rates for each point on the surface, the resulting dis- 
tribution of wate3>-iinping0ment rate over the svirface was obtained. 

The above calculations were made using a tabular form of computation, 
as illustrated in table V. Figure 17 shows the resulting distribu- 
tion of impingement. The total rate of water in5>ingement Mg was 
calculated, using equation (9), by measuring the area under the curve 
shown in figure I7. The value of Mg was cailcxilated to be 2.1 
pounds per hour per foot span. 

Step 2 .— Establish the distribution of heat flow from the 
surface. The distribution of heat flow will depend on the type of 
ice— prevention system to be used. If an electrical system is 
planned, the distribution and intensity, once set, will remain 
xmchanged regardless of variations in flight and meteorologiceLL 
conditions. On the other hand, if the system is to be designed to 
utilize heated air or some other fluid, the distribution of heat flow 
from the surface will depend upon the characteristics of the internal 
flow of the fluid as well as the conditions affecting the external 
heat transfer. If such a system is to be used, calciilation of the 
heat-flow distribution will be rather conq)lex, and it is believed 
that assuming a distribution will provide a good basis for starting 
the ceilculations for design. 

The heat— flow distribution and Intensity used in these calcu- 
lations was estimated, based on dfta presented in references 4 and 
26, to be the heat— flow distribution and intensity of the thermal 
ice— prevention system for the wing of the C— 46 airplane (refer- 
ence 20) at the 8-foot chord station. This distribution, which is 
believed to be representative of a probable thermal system, is shown 
in figure 26. 

Step 3 «— Determine the values of convective heat— transfer 

coefficient. The values of measured convective heat— transfer coeffi- 
cient with the estimated form of the turbulent coefficients shown in 
figure 25 were used. 

Step 4 .— Calculate values of surface temperature in the area of 
water Impingement, using equation (24), such that the values of q at 
einy point are equal to the assumed heat flow. (See fig. 26.) 

3tep5«— From the values of surface temperature ceilc Vila ted in 
step 4, coiBpute the convective heat loss using equation (12). The 
curve of convective heat loss is plotted in figure 26. 
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Step 6 .— Calculate the rate of flow of water aft of the region 
of water interception. This was done "by measuring the area in the 
region of water impingement between the convective heat-flow ciirve and 
the curve denoting heat flow to impinging water (fig. 26) to obtain 
the rate of heat dissipation due to evaporation. The rate of evapora- 
tion was computed from equation (20), and was subtracted from the rate 
of impingement to give the rate of water flow aft of the region of water 
interception. The rate of evaporation was calculated to be 1.6 pounds 
per hoinr* per foot span. Subtracting this value from the rate of water 
striking the surface, 2.1 pounds per hour, foot, the rate of flow of 
water aft, then, is 0.5 pound per hom* per foot span. 

Step 7 «— Determine the wetness fraction and make the proper 
modification to the evaporative factor. Using the curve shown in 
figxrre 21, the wetness fraction for a water flow rate of 0.5 pound 
per hour, foot is 3O percent. It is suggested that the values of 
degree of wetness given in figure 21 be used only to the nearest 
10 percent, since more precise usage is considered to be unwarranted. 

The evaporative factor X was then modified by the 30-P®i*cent 
wetness fraction, so that equation (23) becomes 


X = 1 + 1.12 


( ®3~^°k \ PSL 
\'ts— to]f/ Pi 


For these calculations, the value of Pi was taken as the free— stream 
static pressure, so that 


X = 1 + 1.77 


/ 9s-eok \ 

\ tg— toj^ j 


Step 8 .— Calculate values of surface ten5>erature aft of the area 
of water impingement using equation (25) and the revised value of X 
such that the values of q at any point are equal to the assvuned heat 
flow, (See fig. 26.) 

Step 9 .— From the values of surface temperature calculated in 
step 8, compute the convective heat loss, using equation (12). The 
curve of convective heat loss is shown in figure 26, 

Step 10 .— Extend the calculation of convective heat loss until 
the total rate of evaporation, denoted by the area between the curves 
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of convection and total heat flow (except in the region of in^ingement), 
equals the total rate of water impingement. The rate of evaporation 
was computed from the area "between the two ctirves using equation (20), 
For the case of the 0012 airfoil, the extent of heated region required 
for evaporation of 2,1 pounds per hour, foot span was calcxilated to he 
to 18 percent s/c, which is equivalent to l6,5 percent chord. 

It should he noted that the extent of the heated region can he 
decreased hy increasing the intensity of the total heat— flow distri- 
bution np<^ re— calculating the required extent of heated area. 
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TABIE I.- METEQRQLCXHCAL AKD FLIGHT COMDITIONS FOR VJHICH CCRRESPQHDING DATA 
WERE aBTAIMED FOR THE KACA 0012 ELECTRICALLY HEATED AIRFOIL MODEL 
DURIHG FLIGHT IN NATURAL-ICING CONDITICNS 


Icing 

condi- 

tion 

Flight 

number 

Pacific Standeird 
time 

Liquid- 

water 

content 

(gm/m®) 

Mean- 

effective 

drop 

diameter 

(microns) 

Drop- 

size 

distri- 

bution* 

Free- 

air 

tem- 

pera— 

tvire 

(°p) 

Pres- 

sure 

alti- 

tude 

(ft) 

True 

air- 

speed 

(mph) 

Cloud 

type 

1 

39 

2:03 to 2:05 

0.38 

10 

C 

21 + 

9100 

167 

Stratus 

2 

39 

2:13 to 2:16 

.1+1 

10 

C 

21+ 

8980 

162 

Stratus 

3 

39 

2:19 to 2:22 

.38 

9 

C 

23 

9020 

160 

Stratus 

h 

39 

2:23 to 2:26 

.07 

5 

C 

21+ 

8950 

160 

Stratus 

5 

39 

2:28 to 2:31 

.32 

9 

C 

24 

9010 

157 

Stratus 

6 

^3 

12:27 to 12:29 

.58 

11 

C 

25 

10650 

169 

Stratus 

7 

h9 

1:25 to 1:26 

1,00 

16 

D 

11 

8500 

148 

Cumulus 


* Drop-size distributions defined in reference 17 * 
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TABLE II.- METEOROLOGICAL AHD FLIGHT CONDITIONS FOR WHICH CORRESPONDING DATA 
V/ERE OBTAINED FOB THE NACA 65,2-016 ELECTRICALLY HEATED AIRFOIL MODEL 

DURING FLIGHT IN NATURAL-ICING CONDITIONS 


Icing 

condi- 

tion 

Flight 

niHober 

Pacific Standard 
time 

Liquid- 

water 

content 

(gm/m®) 

Mban- 

effective 

drop 

diameter 

(microns) 

Drop- 

size 

distri- 

bution* 

Free- 

air 

tem- 

pera.- 

ture 

(Of) 

Pres- 

sure 

alti- 

tude 

(ft) 

True 

air- 

speed 

(mph) 

Cloud type 

1 

100 

3:15 to 3:19 

0.26 

13 

D 

19 

10750 

167 

Stratocumulus 

2 

102 

1:30 to 1:35 

.3^ 

18 

E 

26 

11300 

168 

Altostratus 

3 

105 

11:01 to 11:06 

.44 

13 

B 

19 

5100 

157 

Stratocumulus 

k 

105 

11:41 to 11:45 

.60 

13 

A 

19 

5060 

140 

Stratocumulus 

5 

105 

2:32 to 2:36 

.42 

19 

D 

21 

5000 

158 

Stratocumulus 

6 

105 

3:07 to 3:12 

• 3 ^ 

22 

C 

20 

5100 

160 

Stratocumulus 

7 

105 

3:18 to 3:23 

.15 

13 

E 

20 

5300 

162 

Stratocmnulus 

8 

111 

12:32 to 12:36 

.09 

l6 

A 

l6 

12450 

178 

Stratocumulus 

9 

111 

1:08 to 1:12 

.28 

29 

A 

13 

9900 

195 

Stratocumulus 

10 

111 

1:15 to 1:19 

.41 

30 

A 

14 

8900 

160 

Stratoc\anulus 

11 

ll 6 

11:39 to 11:40 

.12 

17 

B 

25 

ll 46 o 

185 

Altocumulus 

12 

116 

11:45 to 11:46 

.12 

16 

B 

26 

11100 

184 

Altociunulus 

13 

116 

11:58 to 11:59 

.11 

15 

B 

26 

11100 

165 

Altocumulus 

14 

117 

12:31 to 12:32 

.35 

13 

B 

24 

11100 

180 

Altostratus 

15 

117 

1:13 to 1:14 

.17 

12 

B 

26 

10540 

165 

Altostratus 


*Drop— size distributions defined in reference 17 . 
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TABLE III,- CCMPAEISQN OF CALCULAOED RATES OF WAOER 

IMPTNCEMEBT ahd evafcration over the ieabing edge 

OF THE NACA 0012 EIECTRICALLY HEAIED AIRFOIL 
MOTEL FOR ICING CONDITIONS OF TABLE I 





Calciilated 

Calculated 




rate of 

rate of 

Icing 

condi- 

tion 

Flight 

nvimiber 

Pacific Stemdard 
time 

water 
In^ilnge— 
ment. Mg 
[ lh/(hr) 
(ft. span)] 

water 
evaporai- 
tlon, Ws 
( lh/(hr) 
(ft. span)] 

1 

39 

2:03 to 2:05 

0.30 

0.17 

2 

39 

2:13 to 2:16 

.35 

. 1^1 

3 

39 

2:19 to 2:22 

.27 

.36 

k 

39 

2:23 to 2:26 

.03 

. 01 + 

5 

39 

2:28 to 2:31 

.19 

.19 

6 

^3 

12:27 to 12:29 

, 6 l 

.62 

7 

h9 

1:25 to 1:26 

l.kl 

. 21 +* 


♦Only leading-edge region heated. Rvtnbewk formed. 
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TABLE 17.- COMPARISON OF CALCULATED RATES OF WATER 
IMPINGEMENT AND ETAPORATTON 07EE THE LEADING 
EDGE OF THE NACA 65,2-016 ELECTRICALLY 
HEATED AIRFOIL MODEL FOR ICING 
CONDITIONS OF TABLE H 


Icing 

condi- 

tion 

Fll^t 

nxnnber 

Pacific Standard 
time 

Calculated 
rate of 
water-drop 
Impinge- 
ment, 

[ lh/(hr) 
(ft. span)] 

Calciilated 
rate of 
water 
evaporSf- 
tlon, Wg 
[ lh/(hr) 
(ft. span)] 

1 

100 

3:15 to 3:19 

0.47 

0.46 

2 

102 

1:30 to 1:35 

.94 

.96 

3 

105 

11:01 to H:06 

CVJ 

lf\ 

• 

.55 

4 

105 

11:41 to 11:45 

.56 

.54 

5 

105 

2:32 to 2:36 

.97 

1.05 

6 

105 

3:07 to 3:12 

1.01 

1.06 

7 

105 

3:18 to 3:23 

.33 

GO 

on 

• 

8 

111 

12:32 to 12:36 

.18 

.20 

9 

111 

1:08 to 1:12 

1.60 

.44* 

10 

111 

1:15 to 1:19 

1.79 

.51* 


*Only leading— edge region heated. Runhack formed. 
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TABLE V.- CALCULATIONS FOR liEIEEMINING DISTRIBUTION OF WATEIl-IMPINCEMENT 
RATE OVER THE SURFACE OF AN NACA 0012 AIRFOIL 


n/m 

(percent) 

a/a.me 

a 

(microns) 

K 

5 

0.23 

2.9 

0.003 

10 

0.44 

5.5 

0.012 

CVJ • • • 

0 CO 0 

0 

65 

1 

027 

3 < 

l.< 

12 , 

o.< 

3 

DO 

.5 

063 

2 ( 

1.1 

18 , 

0 .: 

D 

^5 

.1 

L 33 

H • • 

CVJ CV) 0 

D 

DO 

.0 

553 

c 

2 .' 

34 , 

0.1 
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Flgnire 5*— C— 4-6 test airplane as flown during the winter of 194-6— 4-7 showing the position of the 
airfoil models mounted on the f\iselage. 
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(a) NACA 0012 section mounted on C— 46 fuselage 
for the 1945—46 fll^t tests. 



(b) NACA 65,2-016 section mounted on C— 46 fuselage 
for 1946—47 flight tests. 

Figure 6.— Electrically heated airfoil models used to obtain data 
in natural— icing conditions. 
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3. Plastic— impregnated fabric 

4. Electrical resistance heating stn 

5. Plastic- impregnated fabric 


6. Plastic base 
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9 Electricolly heated guard sections 

to. Electrical leads from resistance heating stops 
II. iftewmg blister on fuselage 
!2. Camera 
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Figure 7.- Cut-eway view of the NACA 65,2-0l6 electrically 
heated airfoil model showing construction details. 
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Figure 9 .— Typical runtack formation obtained on the NACA 0012 airfoil 
model with only leading-edge region heated. 
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Figure 13.— Eunback formation obtained on the electrically heated 
NACA 65,2-016 airfoil model during icing conditions 9 and- 
table II. 
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(a) Icing condition 11, table II. 



(b) Icing condition 13, table II. 
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Figure 15 .— lyplcal data records showing area of drop impingement on the NACA 65 , 2-016 airfoil 
model, and traces of water flow aft of the region of impingement. 
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